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Spinning dust emission: the effect of rotation around a 
non-principal axis 
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ABSTRACT 

We investigate the rotational emission from dust grains that rotate around non- 
principal axes. We argue that in many phases of the interstellar medium, the smallest 
grains, which dominate spinning dust emission, are likely to have their nutation state 
(orientation of principal axes relative to the angular momentum vector) randomized 
during each thermal spike. We recompute the excitation and damping rates associated 
with rotational emission from the grain permanent dipole, grain-plasma interactions, 
infrared photon emission, and collisions. The resulting spinning dust spectra gener- 
ally show a shift toward higher emissivities and peak frequencies relative to previous 
calculations. 
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1 INTRODUCTION 

One of the difficulties in measuring the anisotropies in the 
cosmic microwave background (CMB) is that the interstel- 
lar medium (ISM) also emits microwave radiation through 
several mechanisms. This "foreground" radiation must be 
modeled and subtracted in order to measure the cosmologi- 
cal parameters accurately using the CMB. The standard the- 
ory of ISM microwave emission contains three major emis- 



sion mechanisms (e.g. Tegmark et al.l l2000l : iBennett et al.l 



l2003l : iFraisse et al .' 2009) : synchrotron radiation from rela- 
tivistic electrons spiralling in the Galactic magnetic field; 
free-free radiation from ionized gas; and thermal emission 
from dust grains. These are typically traced by external 
templates: low- freq uency radio maps fo r the synchrotro n 
l|Haslam et al.lll982 ''). Ha for the free- free (Finkbcincr 200^"* 



served as a signal at 14.5 and 32 GHz strongly correlated 
with the diffuse 100 ^ m intensity. It was far too bright to 
be thermal dust and had a fiat spectrum across these bands, 
and low-frequency radio and Ha observations predicted far 
too little synchro t ron o r free- free emission to explain the sig- 
nal. iLeitch et al.l (|l997l ') proposed that the signal originated 
from hot gas at T J5 10® K, w hich could produce f ree-free 
radiation but little Ha; however iDraine fc LazarianI (|l998ah 
showed that this gas would cool rapidly and that keeping it 
hot was energetically unfeasible. Several alternative expla- 
nations have been proposed. Spinning dust emission is due 
to the rotation of small dust grains with permanent elec- 
tric dipole moments. The basic mechanism has been known 
for decades (lEricksonl 1957 Hovle fc Wickramasinghelll97Gl : 



and far-infrared continuum for the dust ((Finkbeiner et al 

Il999|). 

Kogut et al](|l996al lbl) reported a spatial correlation be- 
tween Galactic microwave emission at 31.5, 53 and 90 GHz 
and the thermal infrared continuum from dust. They in- 
terpreted the microwave emission as dust-correlated free- 
free radiation, on top of the Rayleigh- Jeans tail of ther- 
mal emiss ion from dust. Their observations were con- 
firmed by Ide Oliveira-Costa et al.l (|l997), who measured 
the microwave int ensity of the Galaxy at 30 and 40 GHz. 
iLeitch et al.l l| 19971 ') claimed the presence of an " anomalous" 
component of Galactic microwave emission, which they ob- 
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iRouan et al.(l 19921 : iFerrara fc Dettmarlll994h . and was sug- 
gested as an explanation for the anomalous emission by 
iDraine fc LazarianI l|l998br i (hereafter DL98b). Magnetic 
dust emission is due to thermal fluctuations of the magnetic 
dipole moments of g rains including ferromagn etic or fer- 
rimagnetic materials jPraine fc LazarianI 1 19991 ). Hard syn- 
chrotron radiation would be a new synchrotron compo- 
nent from young (recently-accelerated) high energy elec- 
trons, proposed to be strongly correlated with the far- 
infrared emission from dust du e to their common a ssocia- 
tion with recent star formation (jBennett et al.l [20031 ). Both 
the spinning and magnetic dust hypotheses predict an emis- 
sion spectrum that peaks in the microwave (the former due 
to the rotation rates of the smallest grains, and the lat- 
ter due to the gyrofrequency in ferromagnetic materials). 
The hard synchrotron hypothesis is now disfavoured due to 
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the low polarization of the anomalous component observed 
by the Wilkinso n Microwave Anisotropy Probe {WMAP; 
iPage et al.ll2007|). its strong morpholosical correl ation with 
dust maps ( jFinkbeinej |2004| : iDavies et al .' 2006") , and evi- 
dence that the anomalous em ission has a rising spectrum at 
low frequencies (< 20 GHz: | de Oliveira-Costa et alj 1 19991 : 
iFinkbeiner et al.ll2004l : IWatson et aLlbOOSl ). 

A key test to distinguish the various models for anoma- 
lous emission is to construct predicted emission spectra and 
compare them to observations. DL98b computed spinning 
dust spectra for a variety of interstellar environments, ac- 
counting for the main processes that affect grain rotation: 
collisions, grain-plasma interactions, infrared emission, and 
radiation-reaction torque on the grain electric dipole mo- 
ment. Model spinning dust spectra have been used exten- 
sively to test (and in some cases disfavour or rule out) the 
spinning dust hypothesis for the anoma lous emission seen in 
the diffuse high-Galactic latitude ISM (iBennett et al.ll2003l: 
Finkbeinej |2004| : iDobler fc Finkbeinerl l2008l: iGold et all 



20091 ) ■ in the Galactic Plane (e.g. IFinkbeiner et al.ll2004) 



and in dense regions such as molecular clouds ( Finkbeinen 
l2004l : IWatsori et al] l2005l: IC asassus ct al' |2006| . 120081 ) and 
Hll re gions (iDickinson et al., 2007. 2009J, supernova rem- 
nants ( Scaife et al.ll2007l ), planetary nebulae llCasassus et al 



2007 ), and an exte r nal g alaxy fNGC6946: iMurphv et al 



2010l ). iDobler et al] (|2009l ) have even used the anomalous 
emission seen by WMAP in the warm ionized medium 
(WIM; traced by Ha) to test dust models; they observe 
a factor of ~ 3 lower anomalous emission than predicted, 
which they tentatively interpret as due to depletion of the 
smallest dust grains (the polycyclic aromatic hydrocarbons, 
or PAHs) in the WIM. 

Recently, the grain rotation problem has 
been revisited by two theo retical groups. 
lAU-Haimoud. Hirata fc Dickinson] (|2009l . hereafter AHD09) 
constructed a more detailed model of grain rotation, follow- 
ing the angular velocities of grains using a Fokker-Planck 
equation and re-evaluating the rotational excitation and 
damping rates using updated grain properties and a more 
sophisticated model for the grain-plasma interactions. 
They also released a public code, SpDust, to compute 
spinning dust spe ctra for any input ph y sical c onditions and 
grain properties. lYsard fc Verstraetd (|2009l ) presented a 
quantum-mechanical treatment of several of these processes 
and computed the resulting emission spectra. 

The existing theoretical treatments of spinning dust, 
however, still contain a number of simplifying assumptions. 
One of the major uncertainties is the grain size distribution 
and typical dipole moment, however this uncertainty can be 
turn ed into a virtue by using it to constrain dust models 
(e.g. IDobler et al]|2009l ). Additionally, there are uncertain- 
ties in the physics of grain rotation, such as the validity 
of the Fokker-Planck approximation or the assumed prop- 
erties such as the evaporation temperature of departing ad- 
sorbed atoms. Some of these pieces of physics are not readily 
amenable to improvement by theoretical calculations, but 
others are. 

The purpose of this paper is to revisit the assumption 
by DL98b and AHD09 that grains rotate around the axis 
of largest moment of inertia due to internal dissipation pro- 
cesses. We argue in particular that PAHs in the diffuse and 
high UV flux phases are likely to be in a random nuta- 



tion state. This is not a trivial detail: a dust grain rotat- 
ing around a non-principal axis emits at multiple frequen- 
cies, including frequencies well in excess of the instantaneous 
grain angular velocity. The fact that electric dipole emission 
depends on the second derivative of the dipole moment fi 
rather than just /i enhances the importance of these higher 
frequencies|j We show in Section l47l1 that for disc- like grains, 
at fixed angular momentum incorporating a random nuta- 
tion state increases the spinning dust emissivity by roughly 
an order of magnitude. Of course, having a random nutation 
state also modifies the processes that change grain angular 
momenta. We investigate each of the major processes and 
find that the typical grain angular momentum is reduced, 
but still find a factor of 1.6 increase in the peak spinning 
dust emissivity j^, and a factor of 1.3 increase in the peak 
frequency for WIM conditions0 

This paper is organized as follows. Section [2] reviews the 
key parameters of the grain models. Section [3] describes the 
expected rotational state of grains and the formalism used 
in this paper (and in the updated SpDust) for describing 
the grain angular momentum distribution. Section [4] con- 
siders the electric dipole emission from grains rotating in a 
random nutation state. Subsequent sections consider spin- 
up and spin-down processes for the grains, taking account of 
nutation: Section[S]considers grain-plasma interactions; Sec- 
tion[6]considers infrared photon emission; and Section[7]con- 
siders collisions. Predicted spinning dust spectra are shown 
in Section [51 where we also explore the sensitivity to some 
of our assumptions. We conclude in Section (9] 

The physical processes affecting grains in non-uniform 
rotation are very complex, and this paper contains some un- 
avoidably long calculations. The reader interested primarily 
in the results may skip directly from the end of Section [4. II 
to the beginning o f Section [51 

We note that iHoang et alj l|2010l ) have recently com- 
pleted a related analysis in which axisymmetric dust grains 
are followed through a 2-dimensional space of angular ve- 
locities {ljj^^,u}±). Our analyses agree on the basic conclusion 
that allowing grains to rotate around a non-principal axis 
results in an increase in the spinning dust emissivity and an 
increase in the peak frequency. 



2 GRAIN PROPERTIES 

The physical properties of dust grains treated in this pa- 
per are unmodified from the model of AHD09. We briefly 
summarize the key points here, but refer to AHD09 and the 
references therein for details. 



1 lYsard fc Verstraetd l l200gD allowed for an arbitrary nutation 
state, but imposed the assumption that tlie grain dipole moment 
be exactly parallel to the axis of greatest moment of inertia, which 
eliminates three of the four frequencies of emission from an ax- 
isymmetric grain. They also did not re-consider the coUisional 
and plasma excitation and drag coefficients. 

^ For ease of comparison with previous results, our WIM condi- 
tions are those of DL98b: density nj] = 0.1 cm~'^, gas tempera- 
ture T = 8000 K, H ionization fraction ?i(H+)/nH = 0.99, and 
radiation field normalization x = 1- 
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2.1 Size, shape, and charge 



The grain sizes are described by their volume- equivalent ra- 
dius a, defined h y V = |7ra^. The fidu c ial size distribu- 
tion is taken from IWeingartner fc Draind l|2001bl ). We con- 
sider only the carbonaceous grains because they dominate 
the population of the smallest grains (typically we find that 
grains with radii > 12A make no significant contribution to 
the spinning dust emission). 

As in DL98b, the large grains are taken to be spher- 
ical and the smallest grains are taken to be planar, as 
appropriate for PAHs, and assume the transition to take 
place at 02 = 6 A or Nc ~ 100 carbon atoms. For sim- 
plicity we assume the population of planar grains to be 
disc-like (although real PAHs can have much more compli- 
cated geometries), with a disc radius R = {^a^d~^Y^^ ~ 
7 A {a/^kf^, where we used the interlayer separation in 
graphite, d = 3.35 A, to determine the volume-equivalent 
radius. In AHD09, it was found that the treatment of the 
smallest grains as planar was of only minimal importance, 
resulting in ~ 10-20% changes in the emissivity jv near the 
peak of the spectrum. This conclusion was however based 
on the assumption of r otation aroun d the axis of greatest 
moment of inertia (e.g. |Purcelllll979l ). which we argue here 
is not appropriate. Indeed, we find a substantial (typically 
factor of ~ 2) increase in the spinning dust emissivity as a 
consequence of the disc-like geometry of the PAHs. 

The grain charge distribution calculation is unmodi- 
fied from AHD09; it is based on th e treatment of charg- 
ing by electron and ion coll isions jPraine fc SutinI 1 19871 : 
IWeingartner fc Draind [2001al) and photoelectric charging 



band (|Hudgins et al.ll2005l ). or by incomplete hydrogena- 
tion (or superhydrog enation) of the peripheral carbon atoms 



I Weingartner fc Draind 2001al) assuming a standard in- 



terste llar radiation field ( Mezger et al.l Il982l : iMathis et al.l 
Il983t ) re-scaled by an environment-dependent multiplicative 
factor X- 



2.2 Dipole moments 

The grain permanent dipole moment is one of the most un- 
certain properties as it is not constrained by the UV/optical 
abso rption or IR emission data typ i cally used in dust m odel- 
ing IWeingartner fc DraindboOlbl : iLi fc Draindl200ll ). Our 
fiducial model is similar to that of AHD09 in assuming a 
multivariate Gaussian distribution (appropriate for the ran- 
dom summation of many bonds with dipole moments) with 
a root-mean-square value taken from DL98b: the intrinsic 
dipole moment is taken to be {fJ.f)^^^ = /JA^^/^i where A^'at 
is the number of atoms and /? is a normalization factor. 
The fiducial value is 0.38 D; this is highly uncertain, al- 
though we note that it is reasonable for PAHs that lack ex- 
act symmetries, e.g. the N-circumcoronene cation sequence 
(CsaHigN^, a PAH that would have zero dipole moment 
were it not for the single substitution) has a calculated dipole 
moment corresponding to 13 = 0.16-1.1 D depend ing on the 
position of the substitution (|Hudgins et al.lfioosi '). 

For a nonspherical grain it makes sense to consider 
the orientation of the permanent dipole moment relative 
to the axis of greatest moment of inertia; that is, we 
can consider both the in-plane dipole moment ^ip and 
the out-of-plane moment flop- An in-plane dipole moment 
in a PAH could be produced by e.g. nitrogen substitu- 
tion, as suggested to reproduce the location of the 6.2 /im 



(|Le Page et al.l 20031) . An out-of-plan e dipole moment, as as- 
sumed bv lYsard fc Verstraetd (1 20091 ) . requires breaking the 
mirror-plane symmetry of the PAH, e.g. via warping due to 
pentagonal rings as occurs in corannulene, C2oHiol| 

In the absence of a definitive rationale for choosing the 
dipole moment to be in-plane or out-of-plane, we take for 
our fiducial model the isotropic ratio {/iop)-(Mip) = 1:2 (i.e. 
assign the same moment on all three axes). This choice is 
very uncertain, however we find that the resulting spectra 
are only weakly sensitive to it - e.g. for the fiducial WIM 
model, we find only a ~ 12% change in the characteristic 
emitted frequency and a ~ 10% change in the total emitted 
power between the extreme cases of a purely in-plane dipole 
moment and a purely out-of-plane moment. 



3 ROTATION OF A DISC-LIKE GRAIN 

Here we review the formalism to describe the rotation of a 
general axisymmetric grain, and the physics that determines 
the nutation angle distribution. 



3.1 General description 

We focus here on the case of an oblate axisymmetric dust 
grain, i.e. one with principal moments of inertia Ii = I2 < 
I3. For a planar grain, which is a reasonable model for a 
PAH, one has I3 = 2/i|3 In free solid-body rotation, the 
angular momentum L and rotational energy E^ot are con- 
served; this implies that the angle 9 between the grain sym- 
metry axis and the angular momentum vector is also con- 
served. We may choose the z-axis to be along the angular 
momentum vector, so that 8 is one of the Euler angles of 
the grain (see Fig. [1} . The remaining two Euler ang les 0, 
then advance at a rate (e.g. Eqs. 8.46,47 of .Hand fc FinchI 
Il998l ): 



L_ 
h 



and 



-L(Jr^ -/g"^) cose. 



(1) 



(2) 



We note that ip may have either sign, but one always has 

\i>\ < 0- 

The rotational energy is given by 

£rot = ^ - ^(/r' - /a"') cos^ e. (3) 

The quantum mechanical description of the system will 
occasionally be useful (e.g. for counting states) even though 
the calculation of this paper is in the classical regime, as 
DL98b showed that in general L ^ h even for the smallest 



^ We note that searclies for corannulene rotational lines i n tiie 
Red Rectangle have returned null results l lPilleri et al ]|2009l) . but 
this does not rule out larger warped PAHs. 

For warped PAHs, /3//1 is not exactly 2; but it is e.g. 1.93 
for corannulene according to the structural parameters given in 
iHcdbcrg ct al.. (.2000. ). 
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k 

(disc normal) 




Figure 1. The definitions of the Euler angles used in this paper. 
The xyz axes correspond to the inertial frame and the ijk axes 
to the grain frame. The angular momentum vector lies in the z 
direction, the normal to the grain disc lies in the k direction, and 
the grain permanent electric dipole moment lies in the ifc-plane. 

grainsI3 This descrip tion is the sa me as that for a rotating 
oblate molecule (e.g. lKrotolll992l . §3.9): the good quantum 
numbers are the total angular momentum quantum number 
J, with total angular momentum J(J + 1); its projection 
on the grain 3-axis HK; and its projection on the inertial 
frame z-axis hM. These all take integer values, with J 
and |Ar|,|M| ^ J. In the case where the lab frame 2-axis 
is aligned with the conserved angular momentum, we have 
M = J. The nutation angle satisfies 



cos 9 - 



K 



1 



for J > 1 



(4) 



3.2 Rotational configuration 

The rotational state of PAHs undergoing thermal spikes has 
been studied in many previous works, particularly those con- 
cerned with the polarization of the PAH emission bands (e.g. 
lLegej|l988l : ISironi fc Drainell2009l ) . Here we recount the key 
results and explain why we expect PAHs in the diffuse ISM 
phases to generally not rotate around a principal axis of 
inertia. 



3.2.1 Effect of thermal spikes on grain rotation 

The rotational state of an oblate dust grain is generically 
described by both an angular momentum L and the angle 9 
between this angular momentum and the axis of symmetry 
of the grain. For large grains, we expect dissipation to bring 
the grain to the states of minimum rotational energy with at 



^ For the same reason, we neglect issues of nuclear spin statistics 
that can arise at small values of J and K for molecules with 
nontrivial symmetry groups. 



fixed angular momentum, i.e. ^ = or ^ = tt. For the small 
grains that dominate spinning dust emission, however, the 
physics is different because the grain undergoes occasional 
thermal spikes (following absorption of each UV photon) 
followed by cooling into the vibrational ground state. During 
thermal spikes, rapid transfer of energy is expected to occur 
between rotational and vibrational degrees of freedom. This 
results in a probability distribution for 9: 



P{9\L)de oc exp 



-Erot{L, 9) 



g{9\L)d9, 



(5) 



where g{9\L) oc sin 9 is the density of states|f] This leads to a 
maximum entropy distribution P{9\L) oc sm9 in the limit of 
fcTvib 3> -Erot, which holds immediately after a UV photon 
absorption. As the grain cools, Tvib drops. However, as the 
grain cools, the density of vibrational states drops, and at 
some temperature Tfr the vibration-rotation energy transfer 
freezes out. We thus expect that the distribution of 9 after 
a thermal spike freezes out at: 



P{9\L)d9 oc exp 



—Eiot{L, 9) 



kTtr 



g{9\L)d9. 



(6) 



We consider disc-like grains for g ^ 6A (Nc < 100 carbon 
atoms). We note that using the iDraine fc model 
for the vibrational spectrum, the fundamental mode is ex- 
pected to be at hui/k = 70(iVc/100)~^/^ K. The freeze-out 
temperature should be at least a few times greater than 
this, depending on the mode spectrum and strength of an- 
harmonic and vibration-rotation couplings. This is greater 
than the rotational kinetic energy in most of the ISM phases 
(or similar to it for high radiation density environments 
such as PDRs) . Thus we expect that in most environments, 
fcTfr exceeds the rotational energy, and the direction of 
the grain symmetry axis is almost completely isotropized 
[P{9\L) oc sin^] following each thermal spike. 



3. 2. 2 Frequency of thermal spikes 

Given the major effect of thermal spikes on the rotational 
state, it is important to consider how the time between ther- 
mal spikes Tabs compares to the timescale for changes in 
grain angular momentum Trot. The characteristic timescale 
between UV photon absorptions for a grain of volume- 
equivalent radius a is 



Tabs — 



)abs(a; i^)-^Av 
riv 



(7) 



where — X'^'^.isrf is the ambient radiation field and 
Tra'^Qabs is the absorption cross section. 

The characteristic rotational damping (or excitation, in 

steady-state) timescale for such a grain is Trot = -L ^ , 
where L is the characteristic angular momentum of the grain 
and ^ is the rotational damping rate evaluated at L. Eval- 
uating Trot requires an analysis of the rotational dynamics. 
The AHD09 analysis suggests 



^ An easy way to see that the density of states is oc sin 6 
to note that at fixed total angular momentum J, since K 



\J J{J + 1) cos 9, the number of states per unit cos 9 is constant, 
and hence the number of states per unit S is oc sinS. 
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Trot 



F ' 



G 



1/2 



(8) 



(f) = 20, and ip 



-SI cos 



(10) 



where F and G are the normahzed damping and excitation 
rates; and th and Tod are the idealized characteristic damp- 
ing timescales through collisions with hydrogen atoms and 
electric dipole radiation respectively (see AHD09 and the 
next section for their precise definitions; and note that r^d 
is defined for thermally rotating grains, but that the actual 
dipole damping time varies depending on whether rotation 
is sub- or super-thermal) Q 

Since the smallest grains rotate fastest and determine 
the peak of the spinning dust spectrum, we evaluate the 
above timescales at the smallest grain size a — 3.5 A, for the 
idealized interstellar environments defined in DL98b, Table 
1. We show these timescales in Table [Jl for both the case of 
^ = (AHD09) and for isotropized 9 (using the formulae in 
this paper). 

In the diffuse ISM phases (CNM, WNM, WIM), thermal 
spikes occur with a rate at least ~ 4 to 6 times higher than 
the processes that change the grain angular momentum. The 
rate difference is even more pronouced in regions of high ra- 
diation intensity (RN, PDR), where the small grains can 
absorb several hundreds of photons during the time it takes 
to change their angular momentum. Therefore we expect an 
isotropic distribution P{9\L) oc smO in these phases. Note 
that this is not true of regions of lower radiation density 
(DC, MC), where thermal spikes occur every few hundreds 
to thousands of years and Tabs S> Trot. In such cases, other 
processes will dominate the distribution of 9 and the re- 
sult may be in between complete isotropization (as assumed 
here) and perfect rotation around the /a axis (9 = 0; as- 
sumed in DL98b and AHD09). An example of such an in- 
termediate case would be the Ma xwclhan distribution for 9 
l| Jones fc Spitzeill 19671 : II azarian fc RobcrgG.1997. ). 

3.3 Angular momentum distribution 

The previous spinning dust analysis by AHD09 followed the 
Fokker-Planck equation for the probability distribution of 
grains as a function of their angular velocity vector oj. Since 
oj is not conserved for a nonspherical grain, the proper vari- 
able to follow instead is the angular momentum L. However, 
in order to maintain a simple connection to previous work, 
we define the variable: 

This is the angular velocity that the grain would have if it 
were able to dissipate the energy associated with its nuta- 
tion; we note that the magnitude of the actual angular ve- 
locity uj satisfies |a;| jf2|. In this paper, the Fokker-Planck 
equation is constructed in terms of Ct. 

For disc-like grains considered in this paper, with 7i = 
^/s, the rotational rates become: 

^ In Eq. ((Sjl the damping time is typically th/ F when linear drag 
processes dominate. When electric dipole damping dominates, e.g. 
in the WIM, the angular velocity is typically (GTed/Tjj)^'^* times 
the thermal angular velocity Wth = {kT/Is)^/'^ (AHD09). Since 
electric dipole emission torque scales as uj^ instead of oj, the actual 
timescale for electric dipole damping is then rcd(t<j/wth)~^, or 

(THTed/G)l/2. 



These results will be needed repeatedly throughout the pa- 
per. 

3. 3. 1 Form of the Fokker-Planck equation 

Following the treatment of AHD09, we write the general 
Fokker-Planck equation for the equilibrium distribution of 

(9 r 1 1 9^ 

[D'(f2)/„(r^^ 



m\ + 



The Fokker-Planck coefficients are 



[i5"^(n)/.(n)] =0. (11) 



- hm and E ^ (CI) = hm ^ -. (12) 

st^o 5t 5t^o 5t 



Here D denotes the mean drift in CI, and E denotes the 
diffusion coefficient tensor. 

It is important to note that, because of the isotropic 
distribution of the direction of the grain symmetry axis (see 
Section I3.2|l , these coefficients are averaged over the angle 
9. More explicitly. 



lim i^pi(n,6») sinMe 

4t->0 Ot 



(13) 



and similarly for E'^ (Q,). 

We now assume an isotropic medium, which is a good 
approximation so long as we are considering the total inten- 
sity spectrum (small deviations from isotropy would result 
in net polarization, which is not the subject of this paper). 
The drift and diffusion terms can then be decomposed as 



D{n) = D{Q,)en 
and 



(14) 



E(n) = £11(^)60 ® en + E^{Q){1 - en ® en), (15) 

where en is the unit vector in the direction of fJ and 1 is 
the identity matrix. The function D{Cl) then denotes the 
rate of damping of rotation, while E\\{Q,) and E^{Q.) mea- 
sure random excitation of the magnitude and direction of 
the angular momentum vector. AHD09 then show that the 
overall distribution function for Q, satisfies the equation 



where 



(16) 



(17) 



Note that D is simply equal to D if the excitation rates are 
isotropic and independent of Q,. This is true for some of the 
mechanisms described, but plasma excitation in particular 
has nontrivial Q. dependence and here Eq. (|17|) is necessary. 

3.3.2 Excitation and damping coefficients 

The D{Q,) and E\\ {Q,) are sufficient to write the Fokker- 
Planck equation but are nontrivial to interpret and vary 
wildly as a function of grain size. For this reason, DL98b 
introduced dimensionless coefficients F and G that describe 
damping and excitation rates relative to those that one 
would obtain from the ballistic impact of hydrogen atoms 
on an idealized spherical grain. These are, for process X, 
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Table 1. Characteristic timescales for UV photons absorption and rotational damping for idealized interstellar phases. The rotational 
damping time is shown for grains rotating about their axis of greatest inertia ("case 1", as assumed in DL98b, AHD09), and for grains 
which are randomly oriented with respect to their angular momentum ("case 2", the subject of the present work). All values are for the 
smallest grains (a = 3.5 A or Nc = 20). 



Phase 




DC 


MC 


CNM 


WNM 


WIM 


RN 


PDR 


Tabs (sec) 




2.0 X 10" 


2.0 X 10^ 


2.0 X lO'^ 


2.0 X 10^ 


2.0 X lO'' 


2.0 X 10* 


6.6 X 10^ 


Trot (sec) [case 


1] 


1.6 X lO'^ 


9.5 X lO'^ 


1.9 X 10* 


2.8 X 10** 


2.1 X 10* 


7.0 X 10^ 


1.4 X 10*^ 


Trot (sec) [case 


2] 


1.4 X lO'^ 


4.1 X lO'^ 


8.2 X lO'^ 


1.2 X 10** 


9.0 X lO'^ 


6.9 X 10^ 


1.1 X 10*^ 



Fxm^^Dxm (18) 

and 

Gxm = (19) 

where th is the idealized damping timescale (whose precise 
definition is given in AHD09) and T is the gas temperature. 

3.3.3 Fluctuation-dissipation theorem 

In their analysis of spherical grains, DL98b and AHD09 ar- 
gued that processes resulting from interaction with a ther- 
mal bath at temperature Tx (notably plasma drag and ex- 
citation) should obey the fluctuation-dissipation theorem, 
D = l2,0,E\\/2kTx . The equivalent result for excitation and 
damping coefficients is that F = (T/Tx)G. No such result 
can apply here because the randomization of the nutation 
degree of freedom during thermal spikes renders the notion 
of a "thermal" distribution for not internally consistent. 
However, the fluctuation-dissipation theorem's close cousin, 
the principle of detailed balance, can be of some use if one 
computes damping and excitation of the actions |i| ~ hj 
and L cos 6 = HK for individual (J, K) levels, and then av- 
erages the resulting coefficients over K. We will need to use 
this technique to compute the plasma drag on a grain rotat- 
ing around a non-principal axis. 



4 ELECTRIC DIPOLE EMISSION 

Our task in computing the emission spectrum falls into two 
major steps. One is to relate the distribution of rotational 
states /a(n) to the observable emission. The other, harder 
task, is to compute the Fokker-Planck coefficients arising 
from each mechanism. We consider the emission process in 
this section, and then proceed to consider the damping and 
excitation mechanisms in later sections. 

For the case of a grain rotating around a principal axis 
of inertia, the grain merely rotates with constant angular 
velocity and emits monochromatic radiation at frequency 
a;/(27r). Thus in these models (DL98b, AHD09) the emitted 
spectrum from a particular grain is built up from its dipole 
moment and the probability distribution for lj. The non- 
uniform rotation case treated here is more complicated, as 
we will see that four frequencies are emitted. 



4.1 Emission spectrum 

Our first step in the analysis is to consider how the electric 
dipole moment of a grain varies as a function of time. We 



define the i, j, and k vectors to form a grain-fixed basis 
with k along the symmetry axis. Without loss of generality, 
/i may be assumed to be in the plane defined by i and k. 
Then: 



(20) 



where and ^op are the in-plane and out-of-plane compo- 
nents of the dipole moment, respectively. 

We now consider the behavior of the dipole moment rel- 
ative to an inertial coordinate system. We choose the inertial 
z axis to be parallel to the angular momentum; then we de- 
fine the 3x3 orthogonal matrix U with elements U^i — x i, 
and similarly for the other 8 entries. The entries involving i 
and k axe needed here: 



Uxi = COS (j} COS ip — COS 6 sin ip sin 0, 

Uyi — cos '(/> sin (/> -f cos cos S sin ?/;, 

Uzi = sin tp sin 9, 

Uxk = sin sin </>, 

Uyk = — sin 6 cos (f>, and 

Uzk = cos 6. 



(21) 



For our purposes, it is most convenient to express the first 
two of these using the product-to-sum rule: 

1, 



i 0) cos (■(/) — (j)) 



-f(l + cos &) cos(i/; + <;!))] and 
Uyi = i[(l + cose')sin(V'-l-0) 
-|-(1 - cose*) sin(0 - i))]. 



(22) 



The advantage of this formulation is that since tp and 
are constant, we have expressed all required components of 
U as sinusoidal functions of time. Each sinusoidal function 
directly emits a (5-function spectrum at its frequency. One 
can see that the above components of U oscillate with the 
four (angular) frequencies fp, {ipl, (j) + ip, and tj) — tp. From 
Eq. (|20|) we see that the same frequencies are present in /i 
(as observed in inertial coordinates). 

The power emitted by an accelerating dipole is given by 

(23) 

From Eqs. (|20|) and ((22}, we see we may write /i as 

A = |-^A'ip[(l-cos6')(V)-(/))^cos(i/'-<;/>) 
+ (1 + cos e){ip + <j>f cos{tp + 4>)] 
— /iop0^ sin 6 sin 4>?x 
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+ |--Mip[(l + cose){ip + <j>y sin {tp + , 
+ (1 - cos6){,j) - Tpf sin {(I) - jp)] 



y 



-fj-ipip s'mOs'mipz. 



(24) 



We observe that when we average over many cycles of (j) 
and ip, aU terms average to zero except those which can be 
expressed in terms of just sin^ cj or cos^ uj, where u is one of 
(j), ip, 4> — ip of <l) + ip. Each of these terms contributes power 
which is emitted at frequency to. We find that 



• At frequency (p + ip, the emitted power is 
^ ^ M?p(^ + 0)"(l + cosg)^ 

• At frequency (p — ip, the emitted power is 
^ ^ M?p(0-i/-)"(l-cosg)^ 

'j'-i' 

• At frequency (p, the emitted power is 
^ ^ 2^gp(/,Sin^e 

• At frequency \ip\, the emitted power is 



3^ 



(25) 



(26) 



(27) 



(28) 



The overall emitted spectrum from a grain of given an- 
gular momentum L is then obtained by finding the amount 
of power emitted in a range of angular frequencies (u, u+du) 
using both the emitted power for each of the 4 components 
and the probability of that component falling in the range 
{uj,cj + did). Consider for example the cp + ip component. 
Letting u = cp + xp , we can see that u is bounded by: 

^<tj<2^-^, i.e. Q.i^u)^m. (29) 

Within this range, the probability distribution for uj can be 
found using 



UJ = (p + ip — 



L 






h 


(f 


h 



cos 61 = 0,(2 - cos 61). 



(30) 



Since cos 6 is uniformly distributed between —1 and 1 with 
density ^, the probability density for cj is then 



Prob(a;)da; = i(|-f) d^=^^ 



(31) 



and the nutation angle that corresponds to emission at u is 



= arccos 



The overall emission spectrum for the ^p + (p component is 
then Prob(cj) times the power at this component, Eq. (|25p: 
this iflfl 

IiIlo*[1 + {L/h - cj)/{L/h - L/h)f 



{u) = 



12c«(L//i -L//3) 



(33) 



* Note that P(ui) has units of ergs per second per (radian per 
second) per grain. 



12c3n ■ ^ ' 

A similar calculation shows that we obtain the same spec- 
trum from emission at (p — ip; this is to be expected since 
the two components are related by the symmetry 6 -H- tv — O. 
Thus: 



■ 0-'/' 



(35) 



Following the same procedure, we find that the spectrum 
emitted at \ip\ is given by 



filu^[l^oj^/{L/h~L/hf] 

3c3(L//i -L/J3) 
M?pc.^ (1 ~ c.V»^) 

3c3n 



(36) 
(37) 



within the range O^oj^JL/Ji— L//3, i.e. ^ cj ^ fi. 

Finally the (p component is at angular frequency L/Ii — 
2Q, irrespective of 6. As calculated before, the total power 
emitted at this frequency is 4^op(L//i)*/(9c^). Thus the 
emitted spectrum is 



n(^) 



4^4a/ 
9c3 ' 



4Mop^'* 
9c3 



(5 (cj - 2^) . 



(38) 



The total emitted spectrum is then the sum of the 4 
components, Eqs. (|34H38|) . considered only within their re- 
spective range of validity. In the particular case of /i = i/3, 
we see that L/I^ = fl, L/Ii = 2Q, and 



p{oj\n) 



2 

en 



Xn<cj<3n 



2 

Mip 



1 



4 2 

+ gMop5(^^ - 2^2) 



(39) 



where the truth function x is 1 if the subscripted inequality 
holds and otherwise. The total power emitted per grain is 
then 



P _ 20* / 2 ^ 32 2 

^spdust — I 0/iip + "^Mop 



(40) 



This should be compared to 2n*^fp/(3c'') for the case of 
a grain rotating around the fc-axis; for an in-plane dipole 
moment {flop = 0) the emitted power is 5 times higher, 
whereas for an isotropically distributed dipole moment (/ijp : 
flop = 2:1) the emitted power is ~ 10 times higher. 

The emissivity per H atom (units of ergs~^ Hz~^ sr~^ 
per H atom) can then be obtained by integrating over the 
probability distribution for and the grain size distribution: 

1 dngr 



> = 7; I da- 

2 ./ nn da 



la/ 
2 c3 



da 









3 



1 drigr 
riH da 

3 



2 

Mip 



dn 



UJ 



2A'op _ 2 



(41) 
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where u — I-kv and the factor 1/2 comes from muhiplying 
by 2ti (conversion from uj to u) and dividing by 47r (per 
steradian). 

4.2 Radiation-reaction torque 

We also need the torque —Trad radiated by the tumbling 
dipole. This radiation back-reacts on the grain, applying a 
radiation-reaction torque +Ti.ad- The general formula for 
this torque is 

Trad = -^{Ax/i), (42) 

where (...) denotes a time average. Since the rotation of a 
rigid solid body is quasiperiodic, this amounts to first an 
average over (f> and ip; and in our case, also an average over 
cos 6 because of the rapid redistribution of the nutation an- 
gle. Averaging over (j) immediately implies that the x and 
y components of Trad vanish; the ^-component is, after ex- 
tensive but straightforward manipulation of trigonometric 
functions, 



24c3/3/| 



-3(2J? + lih + /|)Mfp " 8/3>op 



+ [-8/iVfp + 4/|(Mfp + 2mop)] cos 261 



-(/i-/3)'(2/i+/3)MfpCos46l 



(43) 



Averaging over nutation angles (by multiplying by | sin 6 
and integrating over < 6' < tt) gives 



_ -2L^(3J3 + 3Jf J3 + 4J|)^?p + 20L^llii% 
The case of interest here is Ji = i/3, for which 



Q.^ f82 2 32 2 



9 



(44) 



(45) 



This compares with —2Q,'^fj,fp/{3c^) for the uniformly rotat- 
ing case. 

We note that Eq. (|45p can also be obtained semiclasi- 
cally by noting that the photons emitted in the and (j)±ip 
frequencies carry z angular momentum of +h per photon, 
while those emitted at the \tp\ frequency carry no z angu- 
lar momentumlfl The ratio of angular momentum radiated 
to energy radiated is thus uj~^ for the 4> and ± -i/i compo- 
nents, and we could have written 



^^ad,2 



{oj)doj. 



(46) 



This argument, combined with Eq. (|39p . confirms Eq. (|45[) . 

Radiation-reaction is implemented in SpDust using the 
electric dipole damping time Ted, defined by DL98b to be the 
radiation-reaction damping time I//(2|Trad,z|) for a grain ro- 
tating at thermal velocity, i.e. with rotational kinetic energy 



^kT, about the axis of greatest inertia. Mathematically: 



dT 



(47) 



^ This can be seen by observing that the dipole components at 
frequencies </) and 4}±ip are rotating in the xy-plane, while that 
at is oscillating along the z-axis. 



Our calculation establishes that the damping time for planar 
axisymmetric grains is 



_ /|c2 /^82 2 32 2 
' 3kT 145^"'"^ 9 



(48) 



5 PLASMA EXCITATION AND DRAG 

Plasma excitation is the random torquing of dust grains 
via their interaction with passing ions; plasma drag is 
the related effect in which a rotating grain spins down 
by transferring its angular momentum to the surround- 
ing plasma. These processes have been previously com- 
puted for uniformly rota ting grains in several papers 
(jAnderson fc Watson|[l993l : DL98b; AHD09). 

We consider first the excitation in terms of the power 
spectrum of the electric field at the position of the grain. 
Then we consider the drag, which is determined using de- 
tailed balance arguments. Finally, we combine this with the 
analysis of ion trajectories by AHD09 to obtain the plasma 
F and G coefficients. 



5.1 Excitation in terms of electric field power 
spectrum 

The (nutation angle dependent) plasma excitation coeffi- 
cient is given by the usual Fokker-Planck rule, 

liE^^{n,e)At = {ALl). (49) 

This may be evaluated to first order in perturbation theory 
by noting that the change in z angular momentum in time 
At is equal to the integral of the dipole torque, 

/■At 

AL, = / {fl^Ey - flyE^)dt, (50) 

"'0 

where E is the ambient electric field. In terms of the rotation 
matrix U, 



Mip / {UxiEy - UyiEx)<lt 



AL, 



+ /iop / {UxkEy — UykEx)dt. 

Jo 

Then the excitation coefficient is: 



(51) 



iiE^^(n,i 



1 / 



(^Mip {U:c,Ey - Uy,Ex)tdt 



nAt 

X / {UxiEy - UyiEx)t'dt' 

Jo 



+Mop / {UxkEy - UykEx)tdt 



r-At 

X / (UxkEy - UykEx)t'dt' 

Jo 

rAt 

+ 2/iip/iop / (UxiEy — UyiEx)tdt 

Jo 

X j\uxkEy~UykEx)t'dt'y (52) 

To simplify this, we need to change variables to r = t—t' and 
define the inertial frame electric field correlation function by 



CE{r) = {Ex{t)Ex{t')) = {Ey{t)Ey{t'))- 



(53) 
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the XX and yy correlation functions are equal by isotropy 
of the plasma, and the mixed components are uncorrelated, 
e.g. {Ey(t)E^{t')) = 0. We further assume that Ce{t) 
at sufficiently long lag times r, which is appropriate for a 
thermalized isotropic plasma. Then if At is long compared 
to the decorrelation time (as required for the Fokker-Planck 
equation to be valid), Eq. (|52p simplifies to 



7|£|l(f7,e) 



AtCe{t) X 

oo 

+ f4p{U,,kit)U^k{t') + Uyk{t)Uyk{t')) 

+ 2^iip^iop((/:rz(t)(7a.ft(t') + Uy^{t)Uyk{t' )) 



(54) 



5.2 Plasma drag 

The evaluation of the plasma drag is more complicated. In 
principle, it is a result of second-order perturbation theory: 
the dipole moment of the grain modifies the trajectories of 
passing ions, and then the modified charge distribution ex- 
erts a torque on the grain with nonzero expectation value F°l 
However, a much simpler method of evaluating the plasma 
drag is to use the principle of detailed balance to relate 
the rate of small changes in L and 6 to the rate of inverse 
changes. This method works in four stages: first, we need 
to obtain the diffusion tensor due to plasma drag in (J, K) 
space (ignoring the thermal spikes); we need to relate the 
damping rate (A J) to the diffusion tensor; and then we need 
to express D in terms of these coefficients. Finally we per- 
form the average over nutation angles (or equivalently, over 
K at fixed J). 



Note that to obtain this equation we used the fact 
that to lowest order, the electric field and the grain 
orientation are independent, so expressions of the 
type (Ey{t)Ey{t')Uj:i{t)Uxi{t')) can be factored into 

{Ey{t)Ey{t')){Uxi{t)U.i{t')). 

We now perform the angle {(p, xp) averages of the corre- 
lation functions of the U matrix elements using their explicit 
expressions from Eq. (|2ip : for example. 



= sm 



'6l(sin(;/)(t) sin (?!)(*')) 

■ e{s\n(t){t)sm[4>{t) - (t)T]) 



1 2 

= 2 ^^'^ 6 cos{(f>T) 



(55) 



These simplifications give 



/3'i?||(n,e) = 



+ 



drCEii 



( 1 — cos f 



's[(0 - i^)'' 



3[(0 + l/')^ 



+/^op / drC£;(r)sin 9cos{<j>r). 



(56) 



A further simplification can be achieved by switch- 
ing from the electric field correlation function to its power 
spectrum, which is easier to compute. The power spectrum 
Psif) at frequency / is related to the correlation function 
via 



Ce{t) cos ujT dr = Pe ( — 



This reduces Eq. (|56|) to a simple sum 



(57) 



IiE^^{n,e) = ^*?p 



(1 — cos t 



+ ^-^±^Pe 



+Mop sin^ ( 



- ^ 
2tv 



(58) 



The excitation coefficient (fl) used in Eq. (|16|) can then 
be obtained by performing the average over nutation angles. 



5.2.1 Diffusion tensor 

The rate of diffusion of a grain in (J, K) space due to plasma 
excitation is described by a 2 x 2 symmetric diffusion matrix. 
We have already computed the component associated with 
J: 



_ d(AJ^) 

E.j., = ^^ = -E^^{n,e). 

There are also the other components: 
d{AJAJ<) 



E.JK ; 
and 

Ekk 



At 



At ' 



(59) 



(60) 



(61) 



We may compute Ejk by methods similar to those used 
to obtain Ejj. The change h/S.K in the projection of the 
angular momentum onto the grain fc-axis is equal to the 
integral of the projection of the torque onto the fc-axif^. 



h/\K = / k- {fix E)dt 



fiipE ■ j dt. 



(62) 



The evaluation of Eq. (|60|) gives 



E, 



2 r 
Mip 

4fi2 



(l + cose)^PK(^ 



-(1 - cos6')^Pb 



2n 



(63) 



We note that Ejk{J,—K) = —Ejk{J,K) since the two 
terms in brackets are switched (recall that if K — >■ —K then 
9 ^ n — 9 and — >• — i/;). 

A similar technique could also be used to compute Ekk', 
however we will not need Ekk in our analysis because this 
does not enter into the equations for Dj. 



^'^ For the same reason, plasma drag can be thought of as the 
result of emission of plasma "waves" whose amplitude is propor- 
tional to fi and he nce whose angular momentum is proportional 
to | |Ragotil2003) . 

In the second equality here, we have used the triple product 



identity a ■ (b X c) = c ■ (a X b). 
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5.2.2 Relation to drag 

The key to computing plasma drag is the principle of de- 
tailed balance. We note that in true thermal equilibrium 
with the plasma, and in the absence of thermal spikes redis- 
tributing K (i.e. we consider only plasma interactions as a 
mechanism of changing J and K), the probability of being 
in the (J, K) rotational level is 



P{J,K) oc (2J-Hl)exp 



-h^J{J+l) 



X exp 



2hkT 



2kT 



(64) 



with the factor 2J +1 representing the M-sublevel degener- 
acy. 

We define T .k' to be the rate at which dust 
grains in the (J, K) quantum state transition to the (J', K') 
state due to plasma excitation. We further define the quan- 
tum number changes A J = J' — J and AK = K' — K , and 
the mean values J, = {J + J')/2 and /sT. = {K + K')/2. The 
principle of detailed balance tells us that 



r,, 



.K'P{J,K) = Vj,,k'^j,kP{J',K') 



(65) 



Assuming (as appropriate for the Fokker-Planck approxima- 
tion) that lAJI, \AK\ < J, we find 



r 



J.K^J'.K' 



2J' + 1 
2J + 1 

X exp 

AJ 



exp 



-ft2[J'(J'-hl)- J(J+1)] 



2kT 



2hkT 



1 + 



J 



hkT 



kT 

We then define the symmetrized rate, 



(66) 



(67) 



defined at either integer or half-integer values of the argu- 
ments depending on whether AJ and AK are even or odd. 
The rate S is symmetric in the sense that Saj,ak{J*, K.^) = 
S-Aj,-A^(J*,A'.). Then 



J.K^J'K' 



Sa.j,ak{J*, Kt) 



1 + 



AJ rjAJ 



2 J 2/ifcT 



h^{lr^ -l3^)KAK 



2kT 



(68) 



Taylor-expanding S and keeping only terms first order in 
A J and AK gives 



Saj 



1 + 



AJ _ r JAJ 

2 J 2/ifcT 



+ 



h^{I^^ -I^^)KAK 



2kT 



+ ^djSAJ,AK{J,K) 
+ ^dKSAJ,AK{J,K). 



(69) 



We may now relate the excitation rates to the sym- 



metrized rate function. Inspection of Eq. (|59|) gives 
Ejj{J,K) = Yl ^-^'^ 

.J,K^J' ,K' 

AJAK 

= ^J^Saj,ak{J,K), (70) 

AJAX 

and similarly for Ejk and Ekk- We may then investigate 
the mean rate of change of J: 

d{AJ) 



At 



^.J,K-*.J' ,K' 



(71) 



Here the contributions from A J, AK and —A J, —AK nearly 
cancel. They differ only due to the presence of first-order 
terms (in AJ, AK) in Eq. (|69)l : these give 



d{AJ) 
dt 



/AJ 
V 2J 



r JAJ 



AJAK 

h^{ir'-I^')KAK \ 
+ 2fcT j^^^'^^ 



+ ^djSAJ,AK{J,K) 
+ ^dKSAJ,AK{J,K) 



1 

2J 



r J 

2hkf 



h^{I^'-I^')K 

i^J.J H r77= JijK 



2kT 



1 dEjj 1 dEjK 

2 dJ 2 dK ' 



(72) 



We thus arrive at the remarkable result that the rate of loss 
of angular momentum due to plasma drag is expressible in 
terms of Ejj and Ejk- Equation (|72|l is the closest that we 
come to a standard fluctuation-dissipation relation. 

5.2.3 Computation of the drag coefficient D 

In order to continue, we recall that we ultimately need the 
function -D(f2), which first requires us to find D{Q,9) and 
its average over nutation angles D{Q). We recall that 

d(A») 



D{n,e) = -- 



dt 



en. 



This is not the same as 
d(Afi) h d{AJ) 



dt 



Is dt 



(73) 



(74) 



where "Af2" is understood as A|i7| and the last equality 
holds in the large J limit. These two quantities are however 
related: 



AD. 



in + Afii 



= n 



+ 2An ■ a. + (Af2)2 - Q. 

1 f2Afl -fl^^ 
8 I" 



Afin {Any 

^ 20.'^ 



^2 



Averaging and taking the time derivative implies 



d{An) _ d{Aa,) 



en + 



E^{n,e) 



(75) 



(76) 



dt dt n " 

where the parallel part of the excitation was cancelled by 
the third term in Ea. (|75|l . Solving for D then gives: 
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At 



(77) 



and averaging over the nutation angle, or equivalently over 
K, gives the coefficient used in the Fokker-Planck equation: 



+ 



(78) 



h\ At / ,^ ' n 

The modified damping coefficient D{fl) of Eq. (|17|l . is then 



h\ dt 



1 dfiii 



n ^ 2 dn ' 



(79) 



or in terms of J, and using the averaged Eq. (|59p i?! 
^{Ejj)k, we have 



73 ^_ / d(AJ) \ {Ejj)k 1 d 

ft \ dt /^^, J 2dJ 



(80) 



Here the averages are taken over the nutation quantum 
number K. It is critical to note here that d/dj is a total 
derivative, i.e., the averaging over K is understood to take 
place before the differentiation. This is because in the def- 
inition (Eq. I17p . D is ultimately constructed out of drift 
and diffusion coefficients D and E that have already been 
Tf-averaged. Thus we cannot replace the last term with the 
average of a partial derivative, {dEjj/dJ)K- 

Equation (|80|l may be simplified by plugging in Eq. (|72p : 



D = 



1 f'^'J \ IT. ^ 

2J + 2WJ 

^ {KE.,k)k~- 



dJ 



1 / dEjK 
'2\ dK 



1 d ,^ , 



(81) 



By symmetry under change of sign of K (i.e. 9 o n — 8), 
we have Ejk{J,-J) = -Ejk{J,J) and Ejj{J,-J) = 
Ejj{J,J). Also, inspection of Ens. (|58|l and (|63|l in the 
= J (6^ = 0) case shows thaO Sj,/( J, J) = Ejk{J,J)- 
Substituting these results into Eq. (|8ip . we find a mass can- 
cellation resulting in 



AhkT 



EjjAK 



/.7 



AJkT 



KEjkAK. 



(85) 



,7 



The detailed balance-derived drag coefficient can be 
generalized for any bath at temperature Tx, and can be 
written in the form: 

s2 1 



D 



2kTx 



E„ 



K h^EjK 
J 



(86) 



This should be compared to the fluctuation-dissipation the- 
orem, D = /357_E||/2fcrx, valid for a grain rotating around 
its axis of greatest inertia. In particular, Eq. (|86p reduces to 
the fluctuation-dissipation theorem in the limit of spherical 
grain {h/h = 1). 

5.3 Computation of G and F coefficients 

We are finally ready to construct formulas for the G and F 
coefficients. It is most convenient to express these in terms 
of the AHD09 excitation coefficients Gp,ahd(^^), which of 
course have already been calculated. Recall that on account 
of the fiuctuation-dissipation theorem we had -Fp,AHD(n) ~ 
Gp^AHoi^)- In all cases, we set Ji — i/s. 

5.3.1 G coefficient 



5.2.4 Nutation angle average 

Our final step in the above analysis is to perform the average 
over nutation states K. We would like to express Eq. (|8ip 
in a form that does not contain any derivatives, since the 
latter tend to be numerically unstable. We begin by making 
the replacement: 



dK, 



(82) 



valid for large values of J (i.e. the classical regime). Each 
of the three derivative-containing terms in Eq. (|8ip then 
simplifies. For example. 



dEjK 
dK 



dEj 



2Jj_j dK 



-dK 



Ejk{J,J)~Ejk{J,~J) 

2J 



(83) 



The last term simplifies as well: 



dJ 



{Ejj)i 



d_ 
dJ 



EjjdK 



2,P 



EjjdK + 



dEj 



2J./_j dJ 



-dK 



+ 



EjjjJ, J) + EjjjJ, -J) 
2J 



(84) 



We recall that the AIID09 excitation coefficient was derived 
by assuming 6 = 0, in which case Eq. (|58p reduces to 

l!E^^{n;e = 0)=filPE(^^y (87) 

The plasma excitation rate, using Eq. (|19p . is then 

Gp..uum = ^,lP. (^) . (88) 

This allows us to express the electric field power spectrum 
Psif) in terms of the AHD09 excitation coefficients: 

PE{f) = ^^Gp.AHD(27r/). 



(89) 



We may now use Eq. (|58p to obtain the plasma ex- 
citation rate for general 9. Recalling that (f> = 20, and 
7/) = —Qcos9, we find 

Gp{n,9) = Il^^i^Gp.AHD[(2 + cose)n] 



^■^ Altliougli we use the plasma drag calculation to prove it, this 
is a general result. If we start at K = J then a small change 
A J - AK = J - K ^ JAe'^/2. Therefore the combination of dif- 
fusion coefficients Ej J — Ejx evaluates to (AJ(AJ — AiC))/At = 
J{AJA9^)/At. But diffusion is a Ai^/^ process so (AJAO^) is 
at least of order At^^^. Hence, taking the limit as At — >■ 0"*", 
Ejj — EjK vanishes. 
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(i + cose*)^ 

+ ^ ^Gp,ahd[(2 - cos^jn] 

2 

+ ^sm"eGp,AHD(2t7). 
Mip 



(90) 



We now average over values of cos 9 between —1 and +1. The 
first two terms give identical contributions, and the last one 
simplifies using (sin^ 6) = | . Thus 



Gp(n) 



an 



3--j Gp.ahdH-^ 



+ 



,AHD (20.). 



(91) 



5.3.2 F coefficient 

A similar technique works for the drag coefficient. We substi- 
tute Eq. (|85l) for D into Eq. (|18|) to obtain an expression for 
Fp{Q,). This in turn depends on the excitation coefficients 
Ejj [from Eqs. ^ and ^] and Ejk [from Eq. ((63|l ]. 
This leads to an expression involving electric field power 
spectra, which we transform into AHD09 coefficients using 
Eq. (|89|l . Converting the integrals over K to integrals over 
cos 6 = K/J and using hj = /sfi, we find a mass cancella- 
tion of prefactors, giving 



2Gp(n) 

1 ^1 

2 



;i-f-cos6l) = 



Gp,ahd[(2-cos 9)n] 



(1 — cos ( 
4 

cos Sd cos S 



-G 



p.AHD 



[(2 + cos6l)n] 



(92) 



This simplifies to 
/■an 



0J\^ do; 

Gp,AHD(c.)^ 



4/i5p 
3m?p 



Gn ,AHD (2fi). 



(93) 



Note that due to the ii'-averaging there is no longer 
a definite relation between F and G. However, since 
Gp,ahd('^) > for all u, we find the inequality 



1 < 



< 3. 



(94) 



The calculation of Gp.ahd(i^) is one of the most time- 
consuming parts of SpDust, but it varies slowly as a func- 
tion of frequency and is only required over a range of a factor 
of 3 in frequency (il < cj < 3^) . Thus we implement it with 
an approximate integrator as described in Appendix lAl 

We show the plasma excitation and drag coefficient for 
a grain containing A^'c ~ 54 carbon atoms (equivalent radius 
a Ri 5A), in WIM conditions (as defined in DL98b) in Fig.[2l 



6 INFRARED EXCITATION AND DAMPING 

Another major spin-up/down mechanism for the smallest 
grains is the emission of infrared photons during thermal 
spikes. Here we consider the excitation and damping due to 
these spikes. 



1.00 



0.10 



Gp, AND (0. e = 0) 

Fp (n, 6 random) 
Gp (n, 6 random) 



0.01 



0.01 



0.10 1.00 

(I3/6 k jy/'Q 



10.00 



Figure 2. Plasma excitation and drag dimensionlcss coefficients, 
for A^c = 54 in the WIM. 



The excitation rate is doubled from the AHD09 treat- 
ment for all grains (spherical or not) due to a previous error 
associated with the emitted photon angular momentum (the 
angular momentum carried away by a photon is \/2 h rather 
than just h). This is the only modification in this paper that 
applies to spherical as well as disc-like grains. 



6.1 Excitation 

Infrared excitation is the random change in angular mo- 
mentum resulting from the fact that each infrared photon 
emitted by the grain carries away some angular momentum. 
In the previous analyses of DL98b and AHD09, it was as- 
sumed that the resulting change in angular momentum had 
variance (Ai^) = (since a photon carries one quantum of 
angular momentum), or on one axis (ALl) = In fact, 
the excitation is twice this as can be seen from either of the 
following arguments: 

• An electric dipole photon has angular momentum quan- 
tum number j = 1, so if it carries off angular momen- 
tum — AL (thereby imparting AL to the grain via back- 
reaction), we have AL^ = jU + ~ 2ft^. 

• The z-component of the angular momentum of the pho- 
ton is ALz = —mh, where m G { — 1,0, 1} is the azimuthal 
quantum number of the emitted photon. Since these three 
possibilities are equally likely for an isotropically oriented 
grain, we see that (AL^) = f 

Both of these arguments show that the infrared excitation 
GiR is twice that reported in AHD09, i.e. 



GiR 



SnhkT 



— av, 



(95) 



where Fv is the spectrum of infrared radiation emitted by 
the grain (in e.g. ergs^^ Hz^^ sr~^). 

The correct e xcitat ion rate was included in 
lYsard fc Verstraetd (|2009l ). however their formalism is 
quite different (e.g. they use J rather than Q. as the 
independent variable) and so the discrepancy appears to 
have not been noticed previously. 
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6.2 Damping 

We next consider the infrared damping rate, which arises 
due to shght preferential emission of positive over negative 
angular momentum photons from a rotating grain. A clas- 
sical model of the effect can be constructed by consider- 
ing oscillators either in the plane of the grain or out of the 
plane. The torque from an isotropic distribution of oscilla- 
tors would correspond to adding | of the in-plane and ^ of 
the out-of-plane result. 

We consider an oscillating dipole p with angular fre- 
quency LJ — 2-KV and amplitude P. In the out-of-plane case, 
this corresponds to a dipole moment: 



p — Psm{Ljt)k 

— P sin(Ljf)(sin6 sin 0, — sin6' cos ( 

The torque on the grain is then 
„ 2 



3c3 



{P X p), 



(96) 



(97) 



where the derivatives are taken in the inertial frame. The 
evaluation of the time average of the z-component of the 
torque is then a straightforward exercise; to lowest order in 
6, we find 



rpOp 

■ — 



2 r-)2 ■ 2 

uj P sm 



(98) 



A similar exercise for an in-plane oscillator gives 



rpip 



(l-f-cos6')^(0-hi/') 



4c3 

+ (1 - cos6')^(0- ?/>) 



If we average these over nutation angles, we get 
and 



(99) 



(100) 



(101) 



In order to calculate damping coefficients, we must sum 
over all the oscillators P that contribute to the infrared emis- 
sion. The total power emitted by this dipole is 

UJ^P^ 



so we make the replacement: 



du —F, = 



with u = 2ni'. The total torque is then 



(102) 



(103) 



(104) 



where F^,^ and F°'^ are the emission spectra contributed by 
the in-plane and out-of-plane modes. The damping coefB- 
cient is 

FiR = ^/ d,.^^^^-^^. (105) 

We can see that there is very little difference between the 
in-plane and out-of-plane mode contributions (a factor of 



|). Assuming the isotropic case where | of the emission is 
in-plane and i is out-of-planj^. we find 



Fm = 



10th r K 
37r/3 Jo ^'^ 



(106) 



This is I times the AHD09 damping coefficient for spherical 
grains. 



7 COLLISIONS 

Collisions of dust grains break down into several cases: the 
grain may be charged or neutral, and the impactor may be 
ionized or neutral. Furthermore, one must consider not just 
the angular momentum imparted by the incoming particle, 
but also how much angular momentum it carries away when 
it evaporates. We denote the damping and excitation rates 
with subscripts i or „ (for ion or neutral impactor) and super- 
scripts or for incoming or evaporative contributions. 
In the 6 = Q case, the grain's geometry is time-stationary in 
the inertial frame, and incoming particles are equally likely 
to impact the grain whether they approach on prograde or 
retrograde trajectories, and hence P^'"' = 0. For the more 
general case, there will be a new Fi'"'' contribution associ- 
ated with the fact that the grain can physically crash into 
passing particles, and this leads to a preference to accrete 
incoming particles on retrograde orbits. 

The general problem is not tractable analytically, so 
we focus first on the case of neutral impactors on neutral 
grains. We then heuristically extend the calculation to the 
more general case. 



7.1 Damping rate: neutral grains, neutral 
impactors 

There are two contributions to the damping rate. The first 
is the evaporative damping, F^"^^ , which arises because 
particles evaporating off the grain surface preferentially 
have positive Lz- The second is a new contribution, F^^^^ 5 
which arises because a grain rotating around an axis other 
than a symmetry axis preferentially collides with incoming 
particles of negative Lz. We consider both in turn. In both 
cases, we assume the grain to be a convex rigid body whose 
surface area element is d5, whose normal vector is n, and 
whose instantaneous angular velocity is uj. 



7.1.1 Evaporation 

We suppose that a particle evaporates from position r on the 
grain surface. This point has a local surface velocity vq — 
u) X r. The local phase-space density of particles evaporating 
from the grain surface is 



f{r, v) = K exp 



in{v — vo)^ 



(107) 



In the case of the PAH bands, it is known that some bands 
correspond to in-plane vibrations and some to out-of-plane; how- 
ever given the small difference between the two cases, we have 
not tracked them separately. 
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for V in the half-space: 

'H = {v eR'^ : {v-vo)-n> 0}. 



(108) 



The normahzation constant K can be found from the re- 
quirement that the rate of coUisions per unit area is equal 
to the rate of evaporation per unit area. The flux of evap- 
orating particles (in particles cm~^ s~^) is obtained by inte- 
grating {v — vo) ■ nf over H, giving; 



n f 2fcTcv 
2 V m 



K = 



1 diVcoii 
S dt ■ 



(109) 



The angular momentum imparted to the grain by an 
individual escaping atom is obtained from Newton's third 
law, AL = — mr x v. For an ensemble of escaping atoms, 
we should write 

d(AL) _ 



dt 



dS I d''v{-rxv)[{v-vo)-n]f{r,v).{im) 



The velocity integral is straightforwardly evaluated using 
the substitution v = vo + u. The result is 



d(AL) 
dt 



= mK j dS 



-{r X V{))2Ti 



I kTcv 
\ m 



[r X u){u ■ 



(111) 



The second integral has an integrand even in u, so its value 
is exactly ^ of the integral extended over all u G R^. The 
resulting integrand is then a quadratic function of the com- 
ponents of u times a Gaussian. Such integrals are easily 
evaluated; in this case, the result is proportional to r x nr i 
But we know that § r x ndS — for any closed surface, so 
the second integral vanishes Therefore, we keep only the 
first integral. Using vo — (jt: x r, we reduce Eq. (|111[) to 



d(AL) 
dt 



m dA''coii 
" S" dt 



j) r X {ijj X r) dS. 



(112) 



In the particular case of a disc-like grain of uniform and 
infinitesimal thickness, the surface average of r x (oj x r) is 
the same as its volume average, which by inspection is the 
angular momentum L divided by the grain mass M . Thus, 



d(AL) 

dt 



m dA'"coii 
M dt 



m dNco 
M dt 



(113) 



(114) 



This does not depend on 9, so the evaporation contribution 
to the damping is not modified from the principal axis case. 
The relation 

for neutral atoms impacting neutral grains remains valid. 

This can be seen from symmetry, since the integral must be 
linear in r and n, and has the symmetry of a pseudovector. 

This is based on the assumption that the evaporation prop- 
erties (e.g. Tcv) are uniform across the grain surface. If this is 
violated, e.g. by catalytic sites for the formation of H2, then there 
can be a systematic torque. This has been previously investigated 
and found to be negligible for the smallest grains (DL98b). 



7.1.2 Incoming particles 

We now require the angular momentum acquired from in- 
coming particles. This is actually very similar to the previous 
calculation, except that the phase space density of incoming 
atoms has zero net velocity. 



( m 



\2-KkT 



3/2 



exp 



' 2kT J 



(116) 



and the relevant region of velocity space is now the com- 
plement of H, i.e. T-C. The angular momentum transfer rate 
is 

d(AL) _ 



dt 



dS 



d^ii {rxv)[-{v-vo)-fi\f{r,v).{117) 



We now Taylor-expand to first order in vq. The zeroeth- 
order term (i.e. for dq = 0) is proportional to § rxh dS — 
and vanishes. There are two possible contributions to the 
first-order term. One arises from the explicit vo in the inte- 
grand. The other arises from the dependence of the integra- 
tion region T-C^ on vo. That is, to first order in vq, 

d{AL) 



dt 



dS 



d V {r X v){vo ■ h) f{r, v) 



+ 



d V (vq ■ n){r x v){—v ■ h) f{r, v) 



,(118) 



where the integration region H'^ is evaluated at vo ~ 0; 
dH'^ is the boundary of H'^; and d^v is the area element 
on the boundary. The boundary at nonzero vo is displaced 
a distance vo ■ n, hence the inclusion of this factor in the 
second term, combined with the area element d^v, is the 
element of volume that is brought inside H'^ due to nonzero 
f oEflThe second term can be seen to vanish because v fi — 
on the boundary dH'^. Therefore this second term may be 
dropped. 

Using the Maxwellian distribution, we may perform the 
velocity integral in the first (surviving) term in Eq. (|118p to 
get 



— — — = -nW d) dS(r X n){vo ■ n). 

Substituting vq = oj x r , we conclude that 



(119) 



d{AL) 
dt 



-n^^^ jdS{rx n)[(r x n) ■ a;]. (120) 

The triple product implies that this zero if oj, r, and n are 
coplanar, which is the case for grains rotating around an axis 
of symmetry. In our case, however, it is nonzero. We note 
that the average value of the dyadic (r x ri)(r x n) over a 
disc is jR'^{ii + jj), where R is the disc radius and ii + jj 
is the projector into the plane of the grain. Therefore, 



d(AL) 
dt 



(121) 



where ojip is the in-plane part of the instantaneous angular 
velocity. It is equal to u>ip — Lip/Ii. Further using Ji — i/3, 
we find: 



d(AL) 
dt 



mkT Tri?" 



The -I- sign for this term arises because for vq ■ n > 0, H'^ 
expands. 
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V riiu 



n I m 



(122) 



where we have used the definition of th and a^^ = |_R^ in 
the last linO 

In order to complete the derivation, we need the mean 
value of Lip over nutation angles and time. We note that 
this mean value must be in the direction of L, and that 

Lip L = L% = sm 9, (123) 

which has mean value Therefore the mean value of Lip 
is |L and we find: 



(in) 



2 n I m 

3 wh V mn 

The total drag coefficient is the sum, 

.(in) , p,(ev) 5 n fl^ 



^ n 



3 Wh V "T-H 



(124) 



(125) 



(in) 



n I m 3 



AS\r' -Uvnf]. (132) 



For a disc, the integral is vri?'* and a^x = |^*, so it follows 
that: 



1 n I m 



2 riH V "T-H 



(133) 



Thus the coUisional excitation rate for incoming particles is 
the same as it is for the case of the grain rotating around a 
principal axis of inertia. 

The calculation for evaporating particles is the same 
except that we replace T — )■ Tcv'- 



2 wh V "^H T 



(134) 



Once again, there is no difference from the case of rotation 
around I3. 



7.2 Excitation rate: neutral grains, neutral 
impactors 

We now consider the stochastic change in angular momen- 
tum due to collisions with incoming particles. This excita- 
tion rate (unlike the damping rate) can be computed at zero 
grain rotation. The impact of a particle with velocity v at 
position r imparts an angular momentum AL = mr x v. 
The stochastic change in angular momentum along the z- 
axis can be written as 

^!^^=m^ jdS j A'v[z-(rxv)f{-v-n)f{r, v).(l2&) 

The triple product can be cyclically permuted to get 
d{AL2) 



At 



(b dS d%{v-qy(v-n)f{r,v), (127) 



where q = z xr. The integration over velocity is a Gaussian 
times a cubic polynomial over a half-space, which evaluates 
to 



1 



V2n V m 



[q + (q • n) ], 



d{AL^) _ m{kT) 



At 



2tv 



n(pdS[q^ + {q-ny]. 



(128) 



(129) 



Now the integrand is a scalar and hence may be evaluated 
in either inertial or grain-fixed coordinates. We choose the 
grain-fixed coordinates. The nutation angle average is then 
equivalent to averaging over the direction of z, which leads 
to the dyadic relation 

(qq) = i(r'l-rr), (130) 
where 1 is the unit dyadic. This implies that 
{g^ + (q-n)^)=r^-i(r.n)^ (131) 
Plugging into Eq. (I129|l and converting to the G-factor gives 

The excitation radius Ocx is the same as in AHD09 when taking 
the limit of infinitesimally thin disks (AHD09 assumed disks with 
a thickness d = 3.35 A). 



7.3 Excitation and damping: charged grain, 
neutral impactor 

The case of a charged grain is different from a neutral grain 
because of the induced dipole attraction between the grain 
and the atom. The interaction potential is given by 



(135) 



We can solve for the critical separation rc at which the in- 
duced dipole attraction overwhelms the thermal energy of 



the gas, i.e. where V{rc 



\kT: 



3kT 



1.5 



0.67A^ 



1/4 



V8000 K 



-1/4 



A. (136) 



In the cases where rc <^ Ocx, the induced dipole attrac- 
tion is a small perturbation and the coefficients J^i'"^ , Gn , 
Fn'^^\ and Gn^'' are unchanged from the case of a neutral 
grain. On the other hand, if rc ^ Ocx, then an incoming par- 
ticle is certain to impact the grain surface if it passes over the 
barrier in the effective potential Vcrr(r) = L^/mr^ -I- V{r), 
irrespective of the details of the asymmetry of the grain (a 
disc-like grain has a quadrupole moment, but at r ^ flex the 
potential is dominated by the monopole charge). Therefore 
in this alternative case, the shape of the grain is irrelevant, 
and we should use the AHD09 values for the coefficients 
Fi'"' and G[r\ 

While Gn"' is the same in both of our limiting cases 
(rc/ocx 3> 1 or ^ 1), Fn'"' is not the same and it is nec- 
essary to interpolate between the two solutions. We must 
have Fi'°V-Pn ^ for rc/acx » 1 and F^'"' f Fn.AHDog 
for rc/flcx ^ 1. A simple heuristic interpolating functiorF^ 
is 



Since Fn,AHD09 (a^) ^'^^ ^'^ ^ '^^^ heuristic prescrip- 
tion for Fn'"' is such that F^'"' in that limit. It is not clear 
whether Fn'"' should tend to zero for rc 3> cicx, since although 
the relative difference in coUisional rates between impactors on 
prograde orbits and retrograde orbits should vanish, the overall 
collision rate increases because of electrostatic focusing. However, 
only the relative contribution Fn'"^/Fn matters so this should not 
be a concern. 



© 2002 RAS, MNRAS OOClTlfSTI 



16 Silsbee, Ali-Haimoud, and Hirata 



^ n 



(in) 



1 



(137) 



In diffuses phases of the ISM, colhsions with neutral im- 
pactors are in general not the dominant rotational damping 
mechanism (see for example Fig. 4 of DL98b). The exact 
shape of the interpolation function is therefore irrelevant 
in these cases. In very specific environments though, for 
example in extremely dense PDRs, and for low enough 
values of the dipole moment, collisions of neutral impactors 
on charged grains may dominate the rotational damping. If 
this is the case, one should be aware that Fn is uncertain 
in the region rc ~ flex and that this uncertainty will 
propagate on the resulting spectrum, as fpcak y^G/Fn 
and jtot oc [G/Fnf. 



For the evaluation of damping and excitation due to 
evaporating particles, there is no such ambiguity over which 
case to take since we found that the excitation is the same 
for both the uniform ^ = rotation (old case) and isotropic 
6 distribution (new case). 

7.4 Excitation and damping: neutral grain, 
charged impactor 

We now consider the case of an ion impacting a neutral 
grain. The analysis of evaporating particles is the same as 
that treated in Sections l7.1l and l7.2l since the ion is assumed 
to recombine on the grain surface and evaporate as a neutral. 

Incoming ions follow a trajectory influenced by the 
dipole moment of the grain, both permanent and induced. 
The characteristic induced dipole energy for a grain with 
radius R is E\d ~ qcR~^ (i-e- the attraction of the ion to the 
mirror charge). For the PAH sequence, R ~ 0.9Nq^ A, so 
we find Eii/k = 1.7 x 10'^ N~^^'^ K. Thus even at iVc ~ 100 
(our largest disc- like grains), the induced dipole energy is 
well above the temperature of the gas even in warm phases 
(WNM, WIM). Therefore to a first approximation we treat 
the probability of an incoming ion striking the grain surface 
as being determined by the dipole interactions rather than 
grain geometry. (Since the polarizability tensor of the grain 
is not isotropic, this is only an approximation.) In this case, 
we are justified in using the AHD09 rates for incoming parti- 
cles. We are thus led to the conclusion that the AHD09 rates 
are applicable to ion impacts on neutral grains, both for the 
incoming coefficients F, G*!'"' and as previously described for 
the evaporation coefficients F, Gj"^^' . 

7.5 Excitation and damping: charged grain, 
charged impactor 

In the case of an ion colliding with a charged grain, the par- 
ticles interact with the Coulomb potential, which has mag- 
nitude 



V{r) = ^ 



(138) 



for single charges (and more for multiple charges). A sim- 
ple calculation then shows that for practical cases with the 
grains that are treated as disc-like (a < 6 A), ISM temper- 
atures in most phases of interest (including warm phases) 



will have |fcT <C qe/a. In this case, the angular momentum 
transferred to the grain by incoming particles is geometry- 
independent: positive grains will receive essentially no im- 
pacts, while negative grains will be impacted by (and ac- 
quire the angular momentum of) any particle that passes 
close enough to the grain. Thus the incoming rates F, Gj'"' 
are left unaffected. The outgoing rates F, G. "^"^ are also un- 
affected: since the particles are neutral when they evaporate 
off the grain, the outgoing rates are as computed in the pre- 
vious section. 



7.6 Summary 

We may now summarize the key differences between our 
investigation and that of AHD09. 

For the case of neutral grains and neutral impactors, we 
have thus found the coefficients: 



Fn = — -Fn,AHD09 and Gn = Gn,AHD09- 

For charged grains and neutral impactors. 



Fi'"' = <! 1 + ^ 



i+( — 

Hex 



and 

Gn ~ G'n.AHDOg. 



(139) 



(140) 



(141) 



The case of ion impacts is left unchanged from AHD09. 



8 RESULTS 

To avoid lengthy repetitions, we will refer to the case where 
grains are spinning around their axis of greatest inertia (as 
treated by DL98b and AIID09) by "case 1" , and to the case 
where the relative orientation of the grain and the angular 
momentum is randomized (as dicussed in the present work) 
by "case 2". 



8.1 Angular momentum distribution 

We saw in Section 14.11 that, at equal angular momentum, 
the total power radiated by a disc-like grain in case 2 was 5 
times (in the case flop = 0) to ~ 10 times (/ifp : /Xop = 2:1) 
higher than the power radiated in case 1. This ratio goes 
even higher as one increases the fj,^p ■ /i?p ratio. However, 
the angular momentum distribution is different in each case, 
and, as P oc L*, the ratio of the total power emitted will 
really be 

p ~ /r2\2 • y^^"^) 

-Teasel \-LJ /easel 

In what follows we show that {I/^)ease2 < {i^)easel. 

First of all, we showed in earlier sections that the damp- 
ing rates are generally higher for grains spinning around a 
non principal axis. This can be understood heuristically as 
follows: for a given angular momentum L, the rotational en- 
ergy Erot{L, 9) as a function of the nutation angle was given 
in Eq. Q. Averaging over angles, we find that 

(S.ot)(L) = £-^-~ ill" - Is") ■ (143) 
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Figure 3. Probability distribution function for the parameter 
O = L/Ig, for a grain of radius a = 5 A, in WIM conditions, with 
fi?^ : /Xqp =2:1, and with dipolc moment per atom /3 = 0.38 
debye. 




1/ (GHz) 

Figure 4. Power radiated by a grain of radius a = 5 A, in WIM 
conditions, with fi?^ : fi'^p =2:1, and with dipole moment per 
atom 13 = 0.38 debye. 



In the case of a disc- like grain, (Is = 2Ji) this is 
{E,^,)iL) = ^£- = ^EML,e = 0). 



(144) 



Therefore, we may expect that, when in contact with a bath 
of a characteristic energy, grains with a randomly oriented 
rotation axis will have an rms angular momentum ~ ^5/3 
times smaller than those rotating around the axis of greatest 
inertia. This is indeed what we found in the case of collisions 
of neutral grains with neutral impactors, or emission of in- 
frared photons, for which we showed that G was unchanged 
but F was increased by a factor of 5/3. We also showed 
that the normalized plasma damping and excitation rates 
satisfied 1 < Fp/Gp < 3. 

More importantly, the characteristic radiation-reaction 
damping time r^d was found to be shorter in case 2. We have 



Tcd{9 random) 
r,d{e = 0) 



2 



15 /^ip ^ 3 f^°P 
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Figure 5. Spinning dust emissivity in WIM environment. 



In the case where radiation-reaction is the dominant rota- 
tional damping mechanism, which is the case for the smallest 
grains in diffuse phases of the ISM, AHD09 showed that the 

Numerically, we have 



rms angular momentum is oc t^^''' 



1/4 



(6^ random) 
(^ = 0) 



1/4 



0.78 
0.66 



Mop = 0, 



2 

Mip 



flip = 2: 1. 



(146) 



From these considerations, we therefore expect that in the 
same environment, the characteristic angular momentum in 
case 2 will be ~0. 66-0. 78 times the one in case 1. 

We show in Fig. [3] the angular momentum distribution 
for a grain of volume equivalent radius a = 5 A, in WIM 
conditions, with fifp : /lop = 2:1, and with dipole moment 
per atom /3 = 0.38 debye. The rms angular momentum in 
case 2 is ~ 0.67 times the one in case 1. 



8.2 Change in emissivity 

At a given angular momentum, the power radiated in case 
2 peaks at a frequency approximately twice higher than the 
power radiated in case 1 (see discussion in Section [4. ip . 

Therefore, and in view of the preceding section, we ex- 
pect that the total power radiated in case 2 will peak at a 
frequency ~ 2 x 0.7 ~ 1.4 times higher and will integrate to 
a total power ~ 10 X (0.7)'' ~ 2 times the power radiated in 
case 1. This is indeed what we find, as can be seen in Fig.U 

The overall spinning dust emissivity follows the same 
trends, as can be seen in Fig. [S]for the WIM, and in Fig. [S] 
for other interstellar environments. 



8.3 Sensitivity to dipole moment orientation 

It is not clear what is the correct assignment for the direc- 
tion of the grain permanent dipole moment relative to the 
principal axes. Here we analyse the effect of the dipole mo- 
ment orientation on the spinning dust spectrum; it appears 
to make only a minor difference in the WIM environment. 

For the smallest grains where radiation-reaction damp- 
ing is most important, we expect (f2^)^/^ oc r^^^ so 



2\l/2 



1/2 



-1/2 



-1/4 



(case 1), 
(case 2). 



(147) 
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Figure 6. Spinning dust spectra for several environmental conditions: dark cloud (DC), molecular cloud (MC), cold neutral medium 
(CNM), warm neutral medium (WNM), reflection nebula (RN) and photodissociation region (PDR). The environments are defined in 
DL98b, table 1. The parameters for the grain size distribution are Ry = 3.1, be = 6 X 10~^ for the diffuse CNM and WNM phases, and 
Ry = 5.5,6c = 3 X 10~^ for the dense DC, MC, RN and PDR. The dashed line is for a spectrum caculated assuming case 1 {9 = 0), 
whereas the solid line is for case 2 (isotropic 8). The shift to higher frequencies and increase in emissivity in case 2 is systematic for all 
environments. We expect that case 2 should be a better approximation in the diffuse and high radiation intensity phases (WIM, CNM, 
WNM, RN, PDR). 



In case 1 the rotation rate is very sensitive to the orientation 
of the dipole moment (only the in-plane component con- 
ributes to the power and the radiation reaction damping). 
Eventually, when the in-plane component becomes small 
enough, radiation-reaction damping becomes subdominant 
and the rms angular momentum will depend only on in- 
teractions with gas or infrared photons. In case 2 however 
the dependence on fifp/fJ,'^ is quite weak, as the out-of-plane 
component contributes to the power and angular momen- 
tum loss. We show the normalized rms angular momenta in 
case 1 and 2 in Fig. [T] Figure [7] also shows an estimate of 
the peak frequency of the emitted power in both cases. 

The total power radiated by one grain, at a given an- 
gular momentum, was given in Eq. 1)40^ for case 2. Taking 
n ~ (D,^)^^'^, and using the above results, we obtain 

(constant (case 1) 

Thus in both cases the total power is very nearly indepen- 
dent of /i?p//i^. In case 1, when fJ,fp/j-L^ — > 0, radiation- 
reaction damping becomes subdominant and the power be- 



comes proportional to fj,fp. These features are shown in 
Fig. El 



9 DISCUSSION 

The purpose of this work was to revisit the assumption of 
DL98b and AHD09 that PAHs rotate about their axis of 
main inertia. The motivation in doing so is that thermal 
spikes following the absorption of UV photons randomize 
the orientation of the grain with respect to the angular mo- 
mentum axis. These absorption events happen frequently 
enough (i.e. on timescale shorter than the timescale for sig- 
nificant changes in the total angular momentum) that we 
expect such a randomization to be effective in most environ- 
ments. Thus we expect the results from this work ( "case 2" ) 
to be a better approximation to diffuse or high-radiation en- 
vironments (CNM, WNM, WIM, PDR, and RN) than those 
from AHD09, which assumed rapid dissipation of the nuta- 
tional energy (S = or "case 1"). However, the new release 
of SpDust allows the user to choose either case; for exam- 
ple, one may wish to explore the range of cases in dark cloud 
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Figure 7. Top panel: normalized rms angular momentum (£7^)^/^ 
as a function of the ratio of in-plane to total dipole moment. 
Bottom panel: estimate of the peak frequency f i/{dP/du)du j P, 
as a function of this ratio. Both are for a dust grain of radius 
a = SA and dipole moment per atom /3 = 0.38 debye, in WIM 
conditions. 



environments where thermal spikes are infrequent, or what 
happens if an as-yet-unidentified dissipational process is ac- 
tive and restores 9 — 0. 

In this work, we showed that, for a given angular mo- 
mentum, the power radiated by a grain in case 2 is ~ 10 
times higher than that radiated by a grain in case 1. This 
is because in case 2, the grain emits at higher frequencies, 
including above the one corresponding to the instantaneous 
angular velocity, as it is not rotating around the axis of 
greatest inertia. 

We evaluated the rotational excitation and damping 
rates in case 2 as a function of grain size and environment 
conditions, and the resulting angular momentum distribu- 
tion. We showed that in a given environment, grains in case 
2 have a lower rms angular momentum than those in case 1, 
by a factor of ~ 0.7. This is due to larger damping rates, in 
particular radiation-reaction damping, in case 2. 

The combination of these results leads to a spinning 
dust spectrum peaking at slightly higher frequencies in case 
2, and a total power approximately twice as large as that 
emitted in case 1. Finally, we showed that the spectrum in 
case 2 is only weakly sensitive to the precise value of the 
M?p : Mop ratio. 
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Figure 8. Total power emitted by a dust grain of radius a = 5 A 
and dipole moment per atom /9 = 0.38 debye, in WIM conditions, 
as a function of the ratio of in-plane to total dipole moment. 



iDobler et al] (|2009t ) found a tension between theoret- 
ical results and microwave observations of the WIM: the 
theory was a factor of ~ 3 larger than the observations, 
and the peak frequency of the spinning dust and its am- 
plitude could not be simultaneously reconciled by changing 
P (the normalization of the dipole moment). By increasing 
the theoretical emissivity and moving its peak to higher fre- 
quencies, our results may worsen this tension. This seems 
likely to strengthen the empirical case for depletion of the 
PAH population in the WIM phase, however there are other 
conceivable explanations for this discrepancy. The random 
walk model for the dipole moment may not apply well to 
the smallest grains (e.g. one could imagine that some of the 
small PAHs have symmetries that guarantee /i. = exactly) , 
or one could imagine extra low-frequency internal degrees of 
freedom which allow the grain to relax to a state where it 
rotates around the axis of greatest moment of inertia. A de- 
tailed exploration of the parameter space (as was done by 
iDobler et al.|[2009l ) is beyond the scope of this paper. 

As a final note, we present some of the remaining issues 
in the treatment of the rotational physics of the smallest 
dust grains: 

• Triaxiality: Many PAHs have triaxial moment of iner- 
tia tensors (e.g. ovalene C32H14, circumpyrene C42H16, and 
their derivatives). This case was not treated in the present 
paper due to its much greater complexity: since the dipole 
moment then depends on elliptic functions of the angle con- 
jugate to the nutation action, a countably infinite number 
of frequencies are emitted. Aside from this aspect, however, 
the underlying formalism in this paper would be applicable: 
the nutation action (rather than hK = 27rL cos 6) would be 
conserved in free rotation and we would average over this ac- 
tion instead of cos6'. The analysis would also break into two 
cases depending on whether the grain lies on the short-axis 
or long-axis side of the separatrix. 

• Impulsive torques: Some of the sources of torque, such 
as ion impacts, impart large but infrequent changes in an- 
gular momentum. This could in principle lead to "rotational 
spikes" analogous to the well-known thermal spikes in the 
grains' internal energy, and would not be treated correctly 
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by the Fokker-Planck equation (which is a diffusive approx- 
imation) 13 

• Ancillary data: We have not fully quantified the un- 
certainties in the ancillary data, such as evaporation tem- 
peratures, the emissivity in the lowest-frequency vibrational 
modes, and the grain charging model (photoelectric and 
electron/ion impact). However, our hope in making the Sp- 
DuST code publicly available is to provide users the flexibil- 
ity to explore deviations from default or fiducial parameters. 
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APPENDIX A: IMPLEMENTATION OF 
PLASMA DRAG INTEGRALS 

Here we describe our implementation of the plasma drag co- 
efficients, Eqs. ((9T|) and (p3)) . in SpDust. These are integrals 
over the Gp,ahd(i^) function, which is itself time-consuming 
to compute. 

If we wish to calculate the integral J w{x)g{x)dx, where 
w{x) is a known weighing function which propeties will be 
discussed later, and the function g{x) is smooth enough on 
the interval of integration that it can be approximated by 
a quadratic polynomial g{x) ^ a + bx + cx^, then we may 
approximate 

J w{x)g{x)dx « A [g{x+) + g{x-)] , (Al) 

where A = ^ J w{x)dx, and {x+,X-} are the solutions of 
the second order system 



x+ + X- = A ^ J xw{x)dx, 
x\ + x^_ = A~'^ J x'^w{x)dx. 



(A2) 



We now turn our attention to the specific cases of Gp (Q) 
and Fp{Q,). With x — uj/^l and the weighing function w{x) = 
(3 — x)^, we get 

Gp{Q,) ^ --^^Gp,AHD(2r2) 
"^Mip 

+ i[Gp,AHD(n+)-f Gp,AHD(n-)], (A3) 

where 

n± = ^ ^ VWS ^ ^ 1.89»}. (A4) 

Similarly, with the weighing function w{x) — x{3~x)^, 
we get 

Fp{^) ~ ^Gp,AHD(2n) (A5) 

■J/^ip 

+ i [Gp,AHD(n+) + Gp,AHD(t^-)] , (A6) 

where 

^ 8 ± x/W S^ ^ {1.18», 2.02!:;}. (A7) 
o 

We have tested the accuracy of the approximate in- 
tegrator and found that the error was less than 1% in 
the regime where Fp,Gp have significant values, i.e. for 
f2 < Qth = \/okT/l3. More precisely, we checked that 

i^^^^")'xminri,^Uo.01 (A8) 



Fp{n) V ' Fp{ntb) 

for grain radii a = 4, 5, 6 A, gas temperatures T = 
50, 500, 5000 K and grain charge Z 
raly for Gp. 
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